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Spontaneous emission of radiation is one of the fundamental mechanisms by which an 
excited quantum system returns to equilibrium. For spins, however, spontaneous emission 
is generally negligible compared to other non-radiative relaxation processes because of the 
weak coupling between the magnetic dipole and the electromagnetic held. In 1946, Purcell 
realised [1] that the spontaneous emission rate can be strongly enhanced by placing the 
quantum system in a resonant cavity — an effect which has since been used extensively to 
control the lifetime of atoms and semiconducting heterostructures coupled to microwave [2] 
or optical [3, 4] cavities, underpinning single-photon sources [5]. Here we report the hrst 
application of these ideas to spins in solids. By coupling donor spins in silicon to a super¬ 
conducting microwave cavity of high quality factor and small mode volume, we reach for the 
hrst time the regime where spontaneous emission constitutes the dominant spin relaxation 
mechanism. The relaxation rate is increased by three orders of magnitude when the spins 
are tuned to the cavity resonance, showing that energy relaxation can be engineered and 
controlled on-demand. Our results provide a novel and general way to initialise spin systems 
into their ground state, with applications in magnetic resonance and quantum information 
processing [6]. They also demonstrate that, contrary to popular belief, the coupling between 
the magnetic dipole of a spin and the electromagnetic held can be enhanced up to the point 
where quantum huctuations have a dramatic ehect on the spin dynamics; as such our work 
represents an important step towards the coherent magnetic coupling of individual spins to 
microwave photons. 

Spin relaxation is the process by which a spin reaches thermal equilibrium by exchanging 
an energy quantum with its environment {Ug being its resonance frequency) for example 
in the form of a photon or a phonon, as shown in Fig. la. Understanding and controlling 
spin relaxation is of essential importance in applications such as spintronics [7] and quantum 
information processing [8] as well as magnetic resonance spectroscopy and imaging [9]. For 
such applications, the spin relaxation time Ti must be sufficiently long to permit coherent 
spin manipulation; however, if Ti is too long it becomes a major bottleneck which limits 
the repetition rate of an experiment, and in turn impacts factors such as the achievable 
sensitivity. Certain types of spins can be actively reset in their ground state by optical [10] or 
electrical [11] means due to their specihc energy level scheme, while methods such as chemical 
doping have been employed to influence spin relaxation times ex-situ [12]. Nevertheless, an 
efficient, general and tuneable initialization method for spin systems is still currently lacking. 
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FIG. 1. Purcell-enhanced spin relaxation and experimental setup, a) By placing a spin 
in a resonant cavity, radiative spin relaxation can be made to dominate over intrinsic processes 
such as phonon-induced relaxation, b) (top) A planar superconducting resonator comprising an 
interdigitated capacitor in parallel with an inductive wire is fabricated on top of Bi-doped ^^Si. A 
static magnetic field Bq is applied parallel to the (x-y) plane of the 50 nm thick aluminium layer 
with a tunable orientation given by 6. (bottom) Magnetic field lines of the microwave excitation 
field Bi generated by the aluminium wire (arrows) are superimposed over the local concentration of 
Bi donors (red), obtained by secondary ion mass spectrometry (SIMS), c) The sample is mounted 
in a copper box thermally anchored at 20 mK and probed by microwave pulses via asymmetric 
antennae coupled with rate ki « K2/5 to the resonator. Microwave pulses at uq of power Pin 
are sent by antenna 1, and the microwave signal leaving via antenna 2 is directed to the input 
of a Josephson Parametric Amplifier (JPA). d) Energy levels of the electron and nuclear spin of 
bismuth donors in silicon (see Suppl. Info.) 
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At first inspection, spontaneous emission would appear an unlikely candidate to influence 
spin relaxation: for example, an electron spin in free space and at a typical frequency of 
ujs/^TT ~ 8 GHz, spontaneously emits a photon at a rate of ~ 10“^^ s“^. However, the 
Purcell effect provides a means to dramatically enhance spontaneous emission, and thus 
gain precise and versatile control over spin relaxation [1], Consider a spin embedded in a 
microwave cavity of quality factor Q and frequency coq. If the cavity damping rate k = ojq/Q 
is greater than the spin-cavity coupling g, the cavity then provides an additional channel for 
spontaneous emission of microwave photons, governed by a so-called Purcell rate [6, 13] 

r - 

^ '"k74 + (52’ 

where 5 = wq — tUs is the spin-cavity detuning (see Fig. la and Suppl. Info.). 

This cavity-enhanced spontaneous emission can be much larger than in free space, and is 
strongest when the spins and cavity are on-resonance (5 = 0), where Pp = /k. Further¬ 

more, the Purcell rate can be modulated by changing the coupling constant or the detuning, 
allowing spin relaxation to be tuned on-demand. The Purcell effect was used to detect 
spontaneous emission of radiofrequency radiation from nuclear spins coupled to a resonant 
circuit [14], but even then the corresponding Purcell rate Pp ~ 10“^® s“^ (or 1 photon emit¬ 
ted every 300 million years) was negligible compared to the intrinsic spin-lattice relaxation 
processes. In order for photon emission to become the dominant spin relaxation mechanism, 
both a large spin-cavity coupling and a low cavity damping rate are needed: in our experi¬ 
ment, this is achieved by combining the microwave conhnement provided by a micron-scale 
resonator with the high quality factors enabled by the use of superconducting circuits. 

The device consists of two planar aluminium lumped-element superconducting resonators 
patterned onto a silicon chip which was enriched in nuclear-spin-free ^®Si and implanted with 
bismuth atoms (see Fig. lb) at a sufficiently low concentration for collective radiation effects 
to be absent. A static magnetic field Bq is applied in the plane of the aluminium resonators, 
at an angle 9 from the resonator inductive wire, tunable in-situ. The device is mounted 
inside a copper box and cooled to 20 mK. Each resonator can be used to perform inductive 
detection of the electron-spin resonance (ESR) signal of the bismuth donors: microwave 
pulses at ujo are applied at the resonator input, generating an oscillating magnetic held Bi 
around the inductive wire which drives the surrounding spins; the quantum huctuations 
of this held, present even when no microwave is applied, are responsible for the Purcell 
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spontaneous emission. Hahn echo pulse sequences [15] are used, resulting in the emission 
of a spin-echo in the detection waveguide, which is amplihed with a sensitivity reaching the 
quantum limit thanks to the use of a Josephson Parametric Ampliher [6] and demodulated, 
yielding the integrated echo signal quadrature Hq. A more detailed setup description can 
be found in [3]. 

Bismuth is a donor in silicon [18] with a nuclear spin / = 9/2. At cryogenic temperatures 
it can bind an electron (with spin S = 1/2) in addition to those shared with the surrounding 
Si lattice. The large hyperhne interaction ■ ~T between the electron and nuclear spin, 
where ^ and 'T are the electron and nuclear spin operators and A/h = 1475 MHz, produces 
a splitting of 7.375 GHz between the ground and excited multiplets at zero magnetic held 
(see Fig. Id for the complete energy diagram [19]). This makes the system ideal for coupling 
to superconducting circuits [20, 21]. At low helds {Bq < 10 mT, compatible with the critical 
held of aluminum) all Amp = ±1 transitions are allowed, mp being the projection of the total 
spin -|- ^) along Bq. Considering only the transitions with largest matrix element, 

resonator A (a;oA/27r = 7.245 GHz, Qa = 3.2 x 10®) crosses the |F, mi?) = |4, —4) y-)- |5, —5) 
transition, whilst resonator B (a;oB/27r = 7.305 GHz, Qb = 1-1 x 10®) crosses the transitions 
|4, —4) y-)- |5, —5), |4, —3) yy |5, —4), and |4, —2) yy |5, —3) (see Figs. 2a and b). 

The echo signal Aq from each resonator as a function of Bq shows resonances at the 
expected magnetic helds, split into two peaks of full-width-half-maximum Au/2tt ~ 2 MHz 
(see Fig. 2a). As explained in [3], this splitting is believed to be the result of strain induced 
in the silicon by the aluminium surface structure, which is non-negligible at the donor 
implant depth of ~ 100 nm. In the following we focus on the lower-frequency peak of the 
|4, —4) yy |5, —5) line which corresponds to spins lying under the wire (see Suppl. Info). 
Over the region occupied by these spins, the Bi held amplitude varies by less than ±2%, 
as evidenced by the well-dehned Rabi oscillations observed when we sweep the power of the 
refocusing pulse Pin at the cavity input (see Fig. 2c), allowing us to determine the input 
power of a TT pulse for a given pulse duration. 

We measure the relaxation time Ti by performing an “inversion-recovery” experiment [22] 
(see schematic, top of Fig. 2a), with the static held Bq aligned along x (6* = 0). An pulse 
hrst inverts the spins whose frequency lies within the resonator bandwidth KA/‘2n = 23 kHz 
(or KB/2n = 68 kHz); note that this constitutes a small subset of the total number of spins 
since ka,b ^ Au. After a varying delay T, a Hahn echo sequence provides a measure of 
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FIG. 2. ESR spectroscopy and Purcell-limited Ti measurement, a) (top) Dominant electron spin 
resonance transitions of the Si:^°®Bi spin system (see Suppl. Info). We employ two resonators, A and B, with 
frequencies 7.246 and 7.305 GHz, which cross up to three spin transitions in the magnetic field range 0-6 mT, 
as seen in the echo-detected magnetic field sweep (bottom panel, vertically offset for clarity). Subsequent 
spin relaxation measurements were made at the magnetic fields indicated by the arrows, corresponding to 
the \F,mp) = |4, —4) o |5,—5) transition for each resonator. The doublet structure of each transition is 
caused by strain exerted by the aluminum film on the donors (see Suppl. Info), b) Cavity linewidths for 
resonators A and B are measured to be 23 and 68 kHz respectively, c) Rabi oscillations are driven by varying 
the cavity input power of the refocusing 5-ps-long tt pulse, d) The inversion-recovery sequence is used to 
measure the spin relaxation time Ti. Spin polarisation is measured with a Hahn echo sequence (50-ps-long 
7r/2 pulse, delay r = 500 ps, and 100-ps-long tt pulse). Rescaled by its value for T Ti, Aq goes from —1 
when the spins are fully inverted to -1-1 at thermal equilibrium. The pulse durations were chosen such that 
only spins within a narrow spectral range were detected, producing a well-defined Purcell-limited Ti (see 
Suppl. Info). Data in this figure were obtained witlvthe static field Bq parallel to the inductor (0 = 0). 




































the longitudinal spin polarization. Fitting the data with decaying exponentials, we extract 
Ti = 0.35 s for resonator A and Ti = 1.0 s for resonator B. 

For a quantitative comparison with the expected Purcell rate, it is necessary to evaluate 
the spin-resonator coupling constant g = jeiF, mF\Sx\F + l,mF — 1) ||(5Bj_||, where 7e/27r ~ 
28GHz/T is the electronic gyromagnetic ratio and (5Bx is the component of the resonator 
held vacuum huctuations orthogonal to Bq (see Suppl. Info of [3]). A numerical estimate 
yields (7o/2vr = 56 ± IHz for the spins located below the resonator inductive wire that 
are probed in our measurements and for 6 = 0. An independent estimate is obtained by 
measuring Rabi oscillations. Their frequency Qr = 2gQ\/^ directly yields go upon knowledge 
of the average intra-cavity photon number n, which can be determined with a ~ 30% 
imprecision from Pin and the measured resonator coupling to the input and output antennae 
(see Suppl. Info). We obtain (7o/2vr = 50±7Hz for resonator A and 58±7Hz for resonator B, 
compatible with the numerical estimate. The corresponding resonant Purcell spontaneous 
emission time is Pp^ = 0.36 ± 0.09s for resonator A and 0.81 ± 0.17s for resonator B, in 
agreement with the experimental values. 

According to Eq. 1, a Purcell-limited Ti should be strongly dependent on the spin-cavity 
detuning. We introduce a magnetic held pulse of duration T between the spin excitation and 
the spin-echo sequence (see Fig. 3a), which results in a temporary detuning 6 of the spins. 
The echo signal amplitude Aq as a function of T yields their energy relaxation time while 
they are detuned by 6. To minimize the inhuence of spin dihusion [22], the spin excitation 
is performed here by a high-power long-duration saturating pulse (see Fig. 3a and Suppl. 
Info) instead of an inversion pulse as in Fig. 2d. As evident in Fig. 3b, we hnd that the 
decay of the echo signal is well ht by a single exponential with a decay time increasing with 
|(5|. The extracted Ti{6) curve (see Fig. 3c) shows a remarkable increase of Ti by up to 
3 orders of magnitude when the spins are detuned away from resonance, until it becomes 
limited by a non-radiative energy decay mechanism with rate F[^^ = 1600 ± 300 s. Given 
the doping concentration in our sample, this rate is consistent with earlier measurements of 
donor spin relaxation times [23], which have been attributed to charge hopping, but could 
also arise here from spatial diffusion of the spin magnetisation away from the resonator mode 
volume. Figure 3c shows that the Ti{6) measurements are in agreement with the expected 
dependence (Fp((5) -f- Fnr)“^, with the only free parameter in this ht being Fnr. 

Having demonstrated the effect of cavity linewidth and detuning on the Purcell rate. 
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FIG. 3. Controlling Purcell relaxation by spin-cavity detuning, a) In-between their 
saturation and subsequent readout, the spins are detuned from the cavity by ^ = applying 

a magnetic field pulse of amplitude Bs, with — 25 GHz/T for this transition and magnetic 

field, b) The decay of spin polarisation is well fit (lines) to exponential decays, with relaxation 
time constants Ti increasing with the detuning, c) Measured Ti as a function of detuning 6. Line 
is a fit with (rp((5) -|- Fnr)”^, yielding Fj;^^ = 1600s. Note that these measurements are taken 
using resonator B and with 0 = vr/d, which results in Ti = 1.7 s at d = 0. The angle 9 changes by 
at most 10% during the magnetic field pulse. 

we finally explore the effect of modulating the spin-cavity coupling constant g. This can 
be achieved by varying the orientation 6 of the static magnetic field Bq in the x-y plane 
(Fig. lb), adjusting the component of the microwave magnetic field (mostly along y under 
the inductive wire) which is orthogonal to Bq. More precisely, g{6) = 'je{F,mF\Sx\F -|- 
l,mF — 1) 6Bfy cos{9y + dBf^ (noting that 6Bix = 0), and we expect 6Biz <C SBiy for 
the spins lying under the wire that are probed in these measurements. This is verified 
















FIG. 4. Dependence of Purcell relaxation on spin-cavity coupling g. a) Rabi oscillations 
(as in Fig. 2c) measured as a function of field orientation 6 (see Fig. lb) are used to extract b) the 
spin-cavity coupling strength g. Its dependence on 9 is fit to the expression in the main text (red 
line); the non-zero value of g{Tr/2) is due to the finite out-of-plane component of the microwave 
magnetic field, c) Inversion-recovery measurements as a function of 9 confirm that the relaxation 
time Ti (see inset) varies as g{9)'^. Solid line is the Purcell formula prediction using the g{9) 
dependence fitted from (b). All data in the figure were taken using resonator B. 

experimentally by measuring the Rabi frequency as a function of 6, as shown in Fig. 4a 
& b, allowing us to extract g{0)/27i = 58 Hz and g{7i/2)/27i = 17 Hz. As expected, we 
measure longer spin relaxation times for increasing values of 9, as shown in Fig. 4c, with the 
relaxation rate scaling as g‘^{9), in agreement with Eq. 1. Overall, the data of Figs. 3 
and 4 demonstrate unambiguously that cavity-enhanced spontaneous emission is by far the 
dominant spin relaxation channel when the spins are resonant with the cavity, since the 
probability for a spin-flip to occur due to emission of a microwave photon in the cavity is 
1/[1 + rNR/rp((5 = 0)] = 0.999, very close to unity. 

At this point it is interesting to reflect on the important fact that the spontaneous emission 
evidenced here is an energy relaxation mechanism which does not require the presence of a 
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macroscopic magnetization to be effective. Under the Pnrcell effect, each spin independently 
relaxes towards thermal eqnilibrinm by microwave photon emission, so that the sample ends 
np in a fnlly polarized state after a time longer than Fp^, regardless of its initial state. This is 
in stark contrast with the well-known phenomenon of radiative damping [24] of a transverse 
magnetization generated by earlier microwave pnlses, which is a coherent collective effect 
nnder which the degree of polarization of a sample cannot increase. We also note that 
had onr device possessed a larger spin concentration, spontaneons relaxation wonld have 
occnrred collectively, manifesting itself as a non-exponential decay of the echo signal on a 
time scale faster than Fp^ [13] and leading to an incomplete thermalization [6, 25]. Snch 
snperradiant or maser emission [26] reqnires the dimensionless parameter C = Ng‘^/[ kAu) 
called cooperativity {N being the total nnmber of spins) to satisfy C* 3> 1 [6, 25, 27], which 
is not the case here becanse of the large inhomogeneons broadening of the spin resonance 
cansed by strain. 

Onr demonstrated ability to modnlate spin relaxation throngh 3 orders of magnitnde 
by changing the applied field by less than 0.1 mT opens np new perspectives for spin- 
based qnantnm information processing: long intrinsic relaxation times which are desirable to 
maximise the spin coherence time can be combined with fast, on-demand initialisation of the 
spin state. We also anticipate Fnrcell relaxation will offer a powerfnl approach to dynamical 
nnclear polarisation [28, 29], for example by tnning the cavity to match an electron-nnclear 
spin flip-flop transition, enhancing the rate of cross-relaxation to pnmp polarisation into the 
desired nnclear spin state [30]. The Fnrcell rate we obtain conld be increased by redncing 
the transverse dimensions of the indnctor wire to yield larger conpling constants, np to 
5 — 10 kHz, bringing the spontaneons emission time below 1 ms (enabling faster repetition 
rates), as well as a higher sensitivity [3], np to the single-spin detection limit. Finally, onr 
measurements constitute the hrst evidence that vacuum fluctuations of the microwave field 
can affect the dynamics of spins, and are thus a step towards the application of circuit 
quantum electrodynamics concepts to individual spins in solids. 
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Supplementary Material: Controlling spin relaxation with a 

cavity 


BISMUTH DONOR SPIN 

The spins used in this experiment, neutral bismuth donors in silicon, have a nuclear spin 
I = 9/2 and an electron spin S = 1/2. The Halmitonian describing the system[l] is 

^/h = B-(7eS®l-7„l®I) + /lS-I, (SI) 

where A/h = 1.45 GHz is the hyperfine coupling between electron and nuclear spins, 
7e/27r = 27.997 GHz/T and 7n/27r = 6.9 MHz/T are the electronic and nuclear gyromagnetic 
ratios. In the limit of a small static magnetic field {Bq < 50mT), the 20 electro-nuclear 
energy states are well approximated by eigenstates of the total angular momentum F = 
S -I- I and its projection, rriF, along the axis of the applied held. These eigenstates can be 
grouped in an F = 4 ground and an F = 5 excited multiplet separated by a frequency of 
(J -|- 1/2) A/h = 7.35 GHz in zero-held, as shown in Figure 1 of the main text. 

For a given weak static held Bq oriented along z, transitions verifying AFAmp = ±1 may 
be probed with an excitation held orientated along x (or y) since their associated matrix 
element (F, \F + 1, mp ±1) = {F,mp\ Sy\F + 1, mp ± 1) has the same magnitude 

as an ideal electronic spin 1/2 transition (m^l Sy \ms') = 0.5. Gharacteristics for the nine 
AFAmp = -|-1 transitions and the nine AFAmp = — 1 transitions are given in Table I for 
Bq = 3mT. Only the ten transitions with a matrix element greater than 0.25 are shown on 
Figure 2a of the main text. Note that the eight hidden transitions with lower matrix element 
are actually degenerate with transitions with stronger matrix element. The transitions 
probed by our resonators are highlighted in red in Table I. 
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Transitions AFAmp = — 1 

|i^, m_p} O l-F + l^mp — 1) 

{F,mF\Sx\F - 1) 

Transitions AFAmp = +1 

T, mp) o l-F + 1, mp + 1) 

{F,mp\ Sx |F+ l,mp + 1) 

Frequency 

(GHz) 

dfjdB 

(GHz/T) 

|4,-4) o |5,-5) 0.474 



7.300 

-25.1 

|4,-3) |5,-4) 0.423 

1 

1 

0.072 

7.317 

-19.2 

|4,-2) |5,-3) 0.372 

1 

CO 

1 

0.125 

7.334 

-13.8 

|4,-1) |5,-2) 0.321 

|4,-2) |5,-1) 

0.176 

7.351 

-8.1 

|4,0) |5,-1) 0.271 

|4,-1) |5,0) 

0.226 

7.368 

-2.5 

|4,l)e^|5,0) 0.221 

|4,0) |5,1) 

0.277 

7.385 

3.1 

|4,2)o|5,l) 0.171 

|4,l)o|5,2) 

0.327 

7.401 

8.7 

|4,3)o|5,2) 0.120 

|4,2) |5,3) 

0.376 

7.418 

14.2 

|4,4)o|5,3) 0.069 

|4,3)o|5,4) 

0.426 

7.435 

19.6 


|4,4) |5,5) 

0.475 

7.452 

25.3 


TABLE I. Relevant Si:Bi transitions and their characteristics for Bq = 3mT. Highlighted in red 


are the transitions accessible to our resonators 


EXPERIMENTAL DETAILS 

Details on the bismuth implanted sample and an extensive description of the setup are 
included in [3]. We present in the following the exact protocols used to acquire the data 
shown in Figure 2 and 3 of the main text. 

Experimental determination of Ti at resonance 

This part aims to explain the inversion recovery sequence presented in Figure 2 of the 
main text. If we rewrite Eq.l of the main text, we can express 

Ti(5) = Ti (0) (1 + 4(57^2) 

with Ti(0) = Since the probed ensemble of spins has a larger linewidth Aca = 

2 MHz than our resonators, the signal emitted during the spin-echo comes from a subset 
of the ensemble of spins, with a frequency spectrum at least as large as the resonator 
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FIG. SI. Excitation pulse bandwidth effect on Ti measurement, (a) Computed pulse bandwidth, 
respectively for a 5(100)-/is vr pulse, in red (blue) incident on a cavity with k/2tt = 23 kHz (green 
dashes). To illustrate the averaging effect of the pulse bandwidth on Ti measurements, the expected 
Purcell Ti curve (black dashes) as a function of spin-cavity detuning is plotted on the right scale, 
with Ti(0) = 0.35s and /t/27r = 23kHz. (b) Ti measurements for two different vr pulse lengths, 
measured on resonance with resonator A. vr = 100/i s (red) yields Ti = 0.35 s, which is in agreement 
with the Purcell rate . vr = 5/rs (red) yields Ti = 0.65 s, a factor 2 away from the accurate value. 


bandwidth. Spins probed at the edges of the bandwidth of the resonator will have longer 
Pnrcell relaxation times: for instance those detnned hy 6 = k have an expected Pnrcell 
relaxation time hve times slower than the Ti time expected at perfect resonance, Fignre Sla. 
The contribntion of those spins with a longer decay time to the signal will resnlt in an 
averaging effect, meaning that the measnred Ti will be erroneonsly longer than predicted. 

In order to snppress this effect, we rednce the bandwidth of the readont seqnence so as to 
collect signal only from spins very close to the resonance. The response fnnction of a pnlse 
of length tp incident on a cavity with bandwidth n at freqnency Uq is expressed as : 

^ yin(t,(^ - c.„V2) ^ X ^ 


tp{uj - CJo) 


tp{uj-uJo) i + 


As shown on Fignre Sla, for the narrowest bandwidth k/27i = 23 kHz of resonator A, 
pnlses of 5ps are heavily hltered by the resonator and have the same bandwidth whereas 
lOOps-long pnlses have a rednced bandwidth of ~ lOkHz. In case of 100 ps-long excitation 
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pulses, the Rabi frequency is such that only spins with |(f|/27r < 5kHz will contribute to 
the signal. This corresponds to a dispersion of only 5% for the expected Purcell relaxation 
times, which is negligible. To illustrate the averaging effect, two inversion recovery curves are 
shown on Figure Sib with readout pulses of 5ps and lOOps. The former yields Ti = 0.65 s, 
which is a factor 2 higher than predicted by the Purcell effect whereas the lattest yields the 
expected value Ti = 0.35 s. 

Thus Figure 2d of the main text shows an inversion recovery sequence that has a readout 
echo sequence with a narrow bandwidth = 100ps, t 7 r /2 = ^ 7 r/ 2 ) to suppress contribution 
from spins with a lower decay rate, and an inversion pulse with large-bandwidth = 5ps) 
in order to maximize the efficiency of the inversion. 


Protocol used to measure spin-cavity detuning dependent relaxation rate 

The goal of this section is to detail the protocol used in Figure 3 of the main text to 
study the Purcell rate dependence on the spin-cavity detuning, S. 

In order to study this dependence, we detune the spins from the cavity by applying a 
magnetic held pulse. This is done in our setup by adding a pulse generator with 50f2 output 
impedance in parallel to the DC supply of one of the Helmholtz coils which have a 1 Hz 
bandwidth. To minimize the effect of transients, buffer times of Is are added after ramping 
the coil up and down. To limit the loss of signal during those buffer times, we use an angle 
6 = 45° and work with resonator B in order to have a longer Ti(0) = 1.68 s. This Ti was 
measured with inversion recovery. All the data presented in this section as well as in Figure 
3 of the main text were done in a separate run. The quality factor of resonator B dropped 
from Qb = 1-07 x 10^ to Qb = 8.9 x 10"^ due to slightly higher losses, yielding the resonator 
bandwidth K/2n = 82 kHz. 

To observe the long relaxation times, such as those measured in Figure 3 of the main 
text, inversion recovery is not an ideal method. Indeed, when the spin linewidth is broader 
(~ x20) than the excitation bandwidth and when the thermalization time is very long, 
one can observe polarization mixing mechanisms [9, 10], spectral and spatial spin diffusion 
being the most relevant to our case, as the system is only constituted from one species. 
If one tries to measure the relaxation from spins that have been detuned by an amount 
6/271 = {us — u:^)/27 t = 3.8 MHz during a lapse of time T with an inversion recovery sequence 
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(Figure S2a), one observes a double exponential relaxation (Figure S2d, green), pointing 
towards the existence of a spin diffusion mechanism. 

Spin diffusion is prevented by suppressing any polarization gradient along the spin line, 
which leads us to use a saturation recovery scheme instead of inversion recovery. The 
simplest saturation recovery scheme (Figure S2b) consists of sending a strong microwave 
tone resulting in the saturation of the line, producing an incoherent mixed state with the 
population evenly shared between excited and ground states. Nevertheless, a relaxation time 
measured with this scheme still yields a double-exponential decay (Figure S2d, orange), with 
time constant similar to the inversion recovery case. This implies that the saturation of the 
line was insufficient. 

To improve the saturation, one can sweep the magnetic held during the saturation pulse 
so as to bring different subsets of the spin line to resonance and realize a full saturation. 
The adopted sweep scheme is shown on Figure S2c. The relaxation curve acquired with such 
a curve is a simple exponential (Figure S2d, blue), indicating the suppression of the spin 
diffusion effect. 

One can further check the quality of the saturation by measuring the polarization across 
the full spin linewidth immediately after saturation. To realize such scans (Figure S2e), we 
apply the relevant saturation pulse at (Uq, then apply a magnetic held pulse Bs = ( 0 ;^—a;o)/7e 
and measure the echo signal Aq{uJs) with a Hahn echo sequence. When no saturation 
pulse is applied, the measured echo signal y4Qo(i^s) is a measure of the full polarization 
— {Szicjs)) = +1 (black curve) and shows the natural spin linewidth. When studying an 
excitation pulse, the polarization of the spins is given by —{Sz{ujs)) = Aq{uJs)/AQQ{ uJa).i 
where Aq{us) is the measured echo signal. Thus —{Sz{ujs)) = —1 ndicates full inversion, 
{Sz{u;s)) = 0 saturation and —{Sz{oJs)) = +1 return to thermal equilibrium. The green, 
orange and blue curves are taken after respectively a vr pulse (a) and a saturation without 
held sweeps (b) and with held sweeps (c). At resonance, one expects a change of Sz from -1 
to -1-1 for a TT pulse and from -1 to 0 for a saturation pulse. Due to the coil transient time, 
all three curves shows a partial relaxation. If the saturation was optimal and no partial 
relaxation was occurring, one should observe 5*2 = 0 for any detuning 6. Among the two 
saturations (b) and (c) studied here, only the last saturation scheme (c) equally saturates 
the line. The basic saturation has a bandwidth 250 kHz and the vr pulse bandwidth is 
similar to the cavity k = 82 kHz. This conhrms that only in scheme (c) can spin dihusion 
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FIG. S2. Spectral spin diffusion. (a,b,c) Ti measurement sequence -when spins are detuned from 
the cavity by applying a magnetic field Bg, providing a detuning 6 = iOg — coq = 2Tr'yesBs, with 
7efr = df/dB^Bo). a) uses a vr = 5ps pulse to realise a so-called inversion recovery sequence, 
b) and c) are saturation recovery sequences: b) uses a 1-s-long strong microwave pulse sent at 
cavity resonance whereas c has in addition a magnetic field scan shown on the bottom part. 
Depicted in orange is the idealized magnetic field profile due to the 1-Hz bandwidth of the coil 
(orange), (d) Ti measurements for sequence a (green), b (orange), c (blue) for 6 = 3.8MHz. 
Fits (black lines): a & b have a double exponential decay whereas c is a simple exponential. We 
attribute this double-exponential decay to spin diffusion, (e) Spectral profiles of excitation pulses 
a (green), b (orange) &: c (blue). The sequence is as follows: send the excitation pulse, detune 
the spins and measure Aq(uJs)- Black line is the reference profile without any excitation pulse, 
yielding reference {Sz{u}s)) = — Aqo{uJs) / Aqq{lJs) ■ When an excitation pulse is sent, one can access 
{Sz{u}s)) = — Aq{uJs)/Aqq{uJs) ■ Note that neither the tt profile or the saturation profile reach either 
the full inversion -|-1 or full saturation 0 at resonance. This is an artefact due to the coil transient 
time. 
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be fully suppressed and yield a simple exponential decay relaxation. This is this last scheme 
that is used to measure the 22 relaxation rates at different detunings 6 of Figure 3 of the 
main text. 

The global £t shown on Figure 3c of the main text is obtained by using equation 
= Fp + Fnr which may be expressed as Ti(0)“^ ^1 + 4 + Fnr to involve 

only experimentally determined parameters. Indeed, k, is precisely determined by measur¬ 
ing the quality factor of the resonator at very low power while Ti(0) is determined by an 
inversion recovery sequence as mentioned above. S has been determined via precise cali¬ 
bration of the coil pulse, thus the only remaining free parameter in the fit is Fnr, yielding 
F^^ = 1600 s. The errors bars come from the accuracy of the relaxation rates fits. 
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